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7. Introduction 


THIs note is concerned with problems analogous to the collapse of a spherical 
bubble in water or the expansion of a spherical cavity due to the pressure of 
an included gas, as in a submarine explosion. In the case where the cavity 
or bubble is situated in liquid extending to infinity, an elegant solution in 
terms of surface harmonics of zero order is given by Lamb (1932). We 
consider here the more general case where the liquid is bounded by concen- 
tric spherical surface where also the inner surface may be collapsible like 
a bubble or be a cavity containing gas, and further the boundaries are 
subject to pressure, and the liquid is subject to the action of external forces. 
This general case can be treated directly by using the Eulerian equations of 
motion as given by Ramsey (1935). But we show that these problems can 
be equally well dealt with by the method given by Lamb. In fact, this note 
arose out of an attempt, during the course of class lectures, to show that all 
the exercises of this type in the text-books of Milne-Thomson (1938, 14 
Examples in Ex. XV, p. 431) and Ramsey (20 examples, pp. 33-39) could be 
solved by the same method. We also show that it can be applied to the 
case of liquid bounded by cylindrical surfaces. 


2. Method of Solution 


In the case of liquid bounded by spherical surfaces, if a, b be the initial 
radii of the inner and outer surfaces and R, R’ their values at time ¢, the 
velocity potential ¢ is given by  =B(t)/r, and the boundary conditions 


at the surfaces, 1iz., — 96/dr = R, R’ for r =R, R’ respectively are satisfied 
by taking B = R*R which is equal to R’?R’ as follows from the equation of 
continuity or from 


R’? — R* =const.; R?R = RR’ (1) 
Thus © 
¢ =R*R/r = R?R'/r (2) 
To this is added the pressure equation 
tp oo)" d¢_ 9 F (t 3) 
po (5) +5 at F i. 
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and (2) and (3) are sufficient for the solution of the problems concerned. 
The additive function F (t) in (3) becomes determinate when the value of the 
pressure at some point of the fluid is given for all values of ¢, and also the 
initial conditions serve to determine the constants of integration that arise 
from the differential equations for R or R’. 


3. Particular Cases 


We now apply the above method to several particular cases which can 
be found as examples in the books by Ramsey and Milne-Thomson. 


(a) Fluid in the form of a spherical shell :— 


Given p = 0 forr = R, p =p, for r = R’: also for t =0, R = R’=0, 
R= 2ne. 


The pressure equation (3) becomes 
Pig 


, ROR’? Ri tae Re 
pitng tian wocnr 








+ F(t) (4) 
Since p = py for r = R’, this determines F (1) and gives 


a R¢R? RERGORR? 2. 
Bote a ang BRD, ROR AIRY page last 





and the condition p =0 for r =R yields 


R“4R’2 12D 2 m2 ma " 
— Pop ——e 1 RA a R hs a R R* = R’R’ ae £ R’? (6) 





Multiplying by 2R’?R’, integrating, and using (2) and the initial data we get 





9 4D? 
— “Po (3 3\ Rieke 13 De 


which gives the value of R’ at any time ¢. In particular when R’ — NC 
R2— 4p, (n'— 1)'8 
30 n— (n8— 1)18 “@B 


(b) Fluid bounded by concentric spherical surfaces and the inner hollow 
filled by zas:~ 
given p = II for r =R’, Q 
obeys Boyle’s law. 





=0; att =0,p =p, forR = ¢,, and the gas 


Using (2), the pr ae 
aye (2), the pressure equation (3) can be written in the form (4) or in 
p_ _ RAR? RER + 2RR? 
p 2 ri = tn =a + F (f t). (4’) 
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The pressure of the gas when the internal radius is R is given by Boyle’s law 
as Pp = prcr®/RE, and since this holds for ail ¢,° the substitution P = Pz for 
r = R in (4) determines F(t) and we get 











pi”, RAR? RR ORR? 3 ayia 

a, ae a a BR ahi Rigen (8) 
The condition p = IJ for r = R’ gives 

He { {RRA RAR HORRA es ay, Pes 

p. 27R4 srhigameags +7 omg 0 tack id Vk Re 
Multiplying by 2R2R and integrating from R =c, to R = cs, we get 

3 ee (RR ae 2c : 
Sa Gad =[7R — R3R2 + PAG Jo R| 
ip R Fe R 1 
=F ~ ctot+ PS log (ce/ey) 
2 ? 


where v = R for R =c, and R =O for R =, and R,’ = value of R’ 
when R = c.. 


Thus we get the value of v from 


ge f Ce as ail Re (9) 


(c) Gravitating fluid between r =a and r =b at t =0:— 
Given p = 0 for+.= R and r = R’ for all ¢; 


R 


Here, aQ= (eer fe 
RK r 
4 R? . 
sor (Be =) + 2mp (R™ — r*) ei 5 (10) 


and the pressure equation (3) becomes 


p R!R2 R°R 
p 


r4 r 








tnil- 


3 
$ 7p eas =) = dap (R’?— r?)+ F (¢) (11) 
Since p = 0 for r = R this gives : 


F ((})}=— (RR + $ R?) + 27 (R’2— R2) 
whence 





— (RR + 3 R*)+ 2mp (R? — R2) 


n= 


Pp RAR? a R?R-+ 2RR? 
p r4 r 


3 
. 7p (7 — ae 2m (R’* — r?) (12) 
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The condition p = 0 for r = R’ now gives 








gan RR! | RRA RE (RRe RD ore RD 
1 R 
— 4p (R?— 5) (13) 


At t=0, R=a, R’=6 and R=0, so that 
2(R ae a® 
o- soe — a(R)y-+ 2p (b?— a2) — 7p (0*- 5) 
which gives. the value of (R), as 


a (R)y= — #79 (2a*— ab— 6) (14) 
From (12) the initial pressure at any point of the fluid is given by 


p a? (R) alg a 
0 F(Ro_ a ()y+ 2p (r*— a) — 479 (5) 





p 
and using (14) this gives after simplification 
Peet” b 
po= #79? (r— a) (b~ r) (4+ 1) (15) 
If we multiply (13) by 2R?R, integrate, and use (1) and (2), we get 
R!R? 





— R3R2 4mp 5 5) 47%p 5 _ 4p "5 5 


ay re (R°R’2— a*b2)=0 





a Re= 4zp R’ {2 (R°— b*)— 5 (R&R’5— a*b?) + 3 (R°— a’)} (15’) 
i154 R3(R’— R) 

(d) Pulsating sphere bounded by infinite liquid:— 

Given R= a+b cos nt, p=0 for r= oo 


The pressure equation is given by 





p—Py_  , R#R?, R?R-+ 2RR? 
=—}—-> + fabeshed 
p r 
and with R= a+ b cos nt, this gives 
sue 
Po— P= aiid & (a+ b cos nt) | b(1— 3 sin? nt)+ a cos nt 
, bsin® at 


“bs “tr (a+b cos nt)? 


At the surface of the sphere r= R, and the pressure there is 
P= Po— p bn® {b (1 — 3 sin? nt)+ a cos nt+- 4b sin? nt} 
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and the least value of this is 
| Proin = Po phn? (a + b). (16) 


The above are typical examples, and others can be solved by the same 
method. 


4. Fluid Bounded by Cylindrical Surfaces 


Let R, R’ be the inner and outer radii of the cylinders. If the velocity 
potential be taken as 


d=B log r 
where r is the distance from the axis, and B a function of the time, then for 


the boundary condition at r= R’, — ep ie R’ for all ¢. Thus we choose 





or 

B=7— RR’, ie:, 
d = —R'R’ log r. (17) 
In place of (1) we have 

R’R’= RR; R?2— R2=a?— B? (18) 
The pressure equation now becomes 

12/2 _ : 

, = ee —(R’R’+ R’’) log r+ F (1) (19) 


if there be no external forces. F (t) is determined by the condition p = J7 
for r = R’ for all ¢, and 


ip 
p 


R’ 


=4 (1) RRR’ RK log * Cy 


If the internal cylinder suddenly collapses we can now find its velocity as also 
the initial pressure at any point of the fluid. 
The condition p= 0 for r=R yields 
i 1 Roe > / Dp? 3/2 R’ 
~=t — ]) R’?— (R’R’+ R’?) 1 : 21 
Salat ) (R’R’ + R) log | (21) 
Multiplying by 2RR and using (18), we have after some simplification 
2) d R’ 
eles =— — (R?R* log — 
ae eae ( : R) 
and integrating 


I R2—— R2R2 log R + const. 
p R 
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If for t= 0, R= 5b, R’ =a the constant is determined, and 
. RAG@eds 
2R 2 Peaks gE ay 
R?R? log (x) 5 (OF RY (22) 
which determines the velocity at the internal surface. 
Equation (20) gives for t= 0, if py= initial pressure at r 


(IT— po) /p=— a(R’), log s 
and from (21) for t= 0, 


] ] as a 
—=—atR’), log 
5 a(R’), log a 


Hence 
log a— log r) 
Py= ll fy eee “eae 
i log a— log bj 
log r— log h 
iT ee eee (23) 
log a—- log b 
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